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Abstract 

Given that it will take quite some time for the Large Hadron Collider (LHC) to 
reach its desired luminosity, it is important to investigate the most favorable scenario 
in which supersymmetry (SUSY) may be discovered at the early runs at the LHC. Our 
aim in this work is to find such a scenario within the gravity mediated SUSY breaking 
(SUGRA) framework and select a class of final states that warrant a discovery at the 
very early runs of the LHC. It turns out that such a situation can be associated with 
a scenario where gluinos are sufficiently light and so are the third generation scalars 
while the first two family scalars are heavy. We find that this can be achieved from 
a high-scale set-up with scalar mass non-universality in the third family and gaugino 
mass non- universality with M3 < Mi,M2. We show that the final state channels which 
are most favorable in such a region of parameter space are 46 -|- E^, Ab + (. + E^ and 
2b + 2i + E^. We also justify our claim by comparing the results with a minimal 
supergravity (mSUGRA) scenario with similar gluino mass. 
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1 Introduction 



With the Large Hadron CoUider (LHC) aheady collected data for the run with a center-of- 
mass energy Ecm=7 TeV, the search for physics beyond the Standard Model has reached a 
new height of excitement. Supersymmetry (SUSY) has been one of the most popular scenar- 
ios in this category, due to its attractive theoretical framework and variety of phenomeno- 
logical features it incorporates. Apart from stabilizing the Higgs sector, it also provides with 
a cold dark matter candidate in form of the lightest SUSY particle (LSP) with i?-parity 
(^(^_]^)3-B+i+25-j cQj^ggj-yation [1,2] . 

Out of different SUSY-breaking models, the most popular one is the minimal super- 
gravity or mSUGRA [3]. Here all the SUSY-breaking parameters are derived from four 
and half parameters, namely, the universal gaugino mass (M1/2), the universal scalar mass 
(mo), the universal trilinear coupling (Aq) all at the GUT scale, tan/3, the ratio of the vac- 
uam expectation values (vev) of the two Higgses and the sign of SUSY-conserving Higgsino 
mass parameter fi. However, proposals have been made to go beyond the universality of 
scalar [4-12] and gaugino masses [13, 14] within the SUGRA framework itself, which are 
strongly motivated from supersymmetric Grand Unified Theory (SUSY-GUT) and low en- 
ergy phenomenological constraints, for example, suppression of the flavor changing neutral 
current (FCNC) or CP-violation. 

With the prospect of LHC running at Ecm= 14 TeV, it is quite likely to shed light upon 
SUSY if it exists at the TeV scale. Whether this is possible, and if at all, then when and 
how, depend strictly on the SUSY spectra that nature has chosen for us. The most likely 
manifestation of SUSY will be the excess in certain final state channels over the standard 
model (SM) background fluctuations. 

Following this, the main trends of the phenomenological studies in search of SUSY in 
context of the LHC have been directed either to a model-based approach where one has a 
specific SUSY model under consideration and then to study the consequences in a collider, 
or to interpret the measured excess in signals to extract the SUSY- parameters such as mass, 
spin etc. using kinematic variables and/or distributions [15]. There have been efforts to study 
the so-called 'LHC inverse problem' [16] as well. While all these are very important and of 
absolute necessity at the same time, perhaps a combined knowledge of all such studies may 
provide us with the most useful hint to unravel SUSY at the upcoming collider experiment. 

In this work, we try to address a relatively simple question - what is the most favorable 
supersymmetric spectra allowed by all the search limits and low energy constraints that can 
leave its imprint during the early run of the LHC and correspondingly what are the final 
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states that are favorable for such a spectra. Our approach here is pragmatic, particularly 
from the experimental point of view, as well as, somewhat model independent to start with. 
We of course, associate our proposed spectra to a SUGRA pattern to justify our claim. 

The answer to the question raised here, may look very simple at the first sight: have the 
minimum possible values of all SUSY particle masses allowed by current bounds. However, 
such a choice may not be consistent with a SUSY breaking scenario. We argue that a spectra 
with light gluinos, light third family scalars, but relatively heavy first two families (that can 
be derived from a high-scale non- universal SUGRA pattern), is one of the most favorable 
cases for the SUSY to be discovered at the early LHC runs. It is even more true when we 
look for the final states in form of 4b + E^, 4b + i + Ejr and 26 + 2£ + Ejr as these arises 
from the decays of the dominant SUSY production processes in such a region of parameter 
space and have little contribution from the SM background. To prescribe these as the golden 
modes for early SUSY discovery, we of course, assume that the machine is tuned properly to 
measure missing energy and tag 6-jets. This is also justified because the main characteristic 
of SUSY-signature is high missing energy which is carried away by the LSP in the i?-parity 
conserving framework. 

We see that, in order to keep the third family squarks much lighter than the first two, 
in our proposed benchmark scenario and to associate them with a SUGRA pattern, we 
require a scalar non-universality at high scale. On the other hand, to keep gluinos light 
and low-lying electroweak gauginos above the LEP limit at the same time, we require a 
gaugino mass non-universality with a heirarchy of M3 < Mi,M2 at the high scale. The 
gaugino mass non- universality of the form M3 < Mi, M2 is achievable within the framework 
of the SUSY-GUT [13,14], while a scalar non-universality of such kind can be motivated 
from string-inspired models with flavor dependent couplings to the modular fields [17,18]. 
To justify our claim, we also compare our results with a mSUGRA parameter point that 
gives similar gluino mass. 

Our paper is organized as follows. In Section 2, we discuss the model under consideration 
and the selected benchmark points chosen for further studies concerning collider signatures. 
In Section 3, we discuss the final states that we look for, the details of the collider simulation 
strategy adopted here and the numerical results obtained from this analysis. We conclude 
in Section 4. 
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2 Model, Formalism and Benchmark Points 



In this section we first advocate the favorable SUSY spectra under consideration. We also 
motivate the final state signals in which such a parameter space can be better observed. 
Then we show that such a spectra can be associated to a high-scale non-universal gaugino 
and scalar mass set-up in a GUT- based SUGRA framework. At the end of this section, we 
discuss the benchmark points chosen for studying the collider signals. 

It is easy to appreciate that a scenario which claims to be favorable for discovery should 
have very light gluinos simply because it has a very high production rate at the LHC (gluinos 
are colour octets and the main production through gluons gets enhanced because of high 
gluon fiux at this energy regime). What about the squarks ? If these are light, they also 
compete in production rate with the gluinos, given the fact that they also have strong inter- 
action and comes with three families. In a generic sense, then light squarks are also favorable 
to add to the SUSY-final states. We need to remember also that the gluino dominantly de- 
cays through the on-shell squark-quark giving rise to more likely a jet-rich final state if first 
two generation squarks participate in the decay chain. The jetty final states have a large 
SM backgrounds coming from QCD processes. Associated leptons may arise in such cases, 
but, mostly from the decays of the electroweak gauginos that appear in the cascades in such 
situations. Increasing the possibility of final state leptons is better to see SUSY signals 
over the background fiuctuations in a hadron collider machine. Although the situation with 
gluinos and squarks all being light, can in principle, warrant a discovery at the early run of 
LHC by suitable choice of event selection criteria, the question is can we think of something 
more distinctive? Actually we can exploit the fact that heavy quarks, namely the top and 
the bottom have the property of decaying to leptons or get identified by b-tagging itself. 
So, if the squarks belonging to the the first two generations be very heavy, then the gluinos 
will decay dominantly through sbottom-bottom or stop-top. While, with a similar sbottom 
and stop mass, the gluino is more likely to decay through sbottom-bottom simply because 
of the phase space. Both cases give rise to bottom quark rich final states. Then, if we look 
for the final states with multiple bs, namely, 46 + E^, 4b + i + E/r or 26 + 2£ + E/r, they 
capture the gluinos decaying through sbottom-bottom or stop-top, with the sbottom decay- 
ing through bottom-neutralino or a top-chargino and the stop decaying to top-neutralino or 
bottom-chargino. In each case, we expect bottom reach final state. Also contribution from 
stop and sbottom pair productions contribute to such kind of final states. 

We discuss the final states with the cuts in details in Section 3. 

From the model point of view, it also turns out that both FCNC and CP-violation 
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constraints may be best tackled if one assumes the first two generations of scalars to be 
multi-TeV and (quasi-) degenerate in masses. ^ This so-called 'inverted hierarchy' is favored 
from low-energy constraints as well. 

A scalar mass hierarchy of this sort can be achieved from a high-scale non- universality in 
first two family scalars with the third one. More specifically, if the squark masses get gen- 
erated from two different uncorrelated mass parameters namely, rriQ^ and mg^'^ respectively 
for the third generation and first two generations at the high scale, then choosing a high 
value of mo^'^ and a small value of mQ^, given a particular value of high scale gaugino mass 
parameter can yield a spectra of the pattern discussed above. Although, it is a phenomeno- 
logical framework that we discuss here, it can be motivated from string inspired models with 
fiavor dependent couplings to the modulaii fields [17, 18]. We would like to note that such a 
set up has been discussed early in some articles for studying the collider signature in context 
of the LHC [4,12], but as mentioned early, not specifically in this context. 

In order to have the gluinos as light as possible and at the same time the low-lying 
charginos and neutralinos above the threshold, we require a non-universality in gaugino 
masses to accommodate it in a SUGRA pattern. The hierarchy required at the high-scale 
is M3 < Mi,M2. This can be generated within the framework of SUSY-GUT with an 
underlying SU (5) or 50(10) gauge symmetry [13,14]. Often one can incorporate a dimension 
five operator in the non-trivial extension of the gauge kinetic function /a^($"') in terms of 
the non-singlet chiral superfields of the form 

ReM<j>)F^,F^^'' = '^Tr{F,,<l>'' rn (1) 

where $^ belongs to the symmetric product of the adjoint representation of the underlying 
gauge group as 



SU{5) : (24 X 2A),ymm = 1 + 24 + 75 + 200 (2) 
SO(IO) : (45 X A5)symm = 1 + 54 + 210 + 770 

Gaugino masses become non-universal if the Higgses responsible for the GUT-breaking, be- 
long to the possible non-singlet representations, unlike the minimal supergravity (mSUGRA) 

^We remind the reader that satisfying constraints imposed by electric dipole moments of electron and 
neutron would require very large scalar masses if we like to have finite values for the CP-violating SUSY 
phases. 
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Representation 


Ms: M2: Ml at Mgut 


75 of SU{5) 


l:3:(-5) 


200 of SU{5) 


1:2:10 


770 of 50(10): H SU{A) x SU{2) x SU{2) 


1:(2.5):(1.9) 



Table 1: Non-universal gaugino mass ratios for different non-singlet representations belong- 
ing to SU{5) or 5*0(10) GUT-group that gives rise to the hierarchy of < Mi,M2 at the 
GUT scale. 

framework [13,14]. Interestingly, the representations 75 and 200 belonging to SU{5) or 770 ^ 
of 50(10) break the GUT group to the SM, we obtain the required hierarchy of M3 < Mi, M2 
at the GUT scale. We tabulate the non-universal gaugino mass ratios for these representa- 
tions in Table 1. A linear combination of these non-singlet representations with the singlet 
one can yield the exact hierarchy that we use in the benchmark points. We would also 
like to mention that such a gaugino mass hierarchy is also supported from the dark- matter 
consideration, as the immediate effect of having a smaller M3 yields a smaller /i after the 
RGE running, yielding a more Higgsino like lightest neutralino, which has much efficient 
annihilation rate to yield a consistent cold dark matter relic density. 

Now we discuss the benchmark points chosen for the study of collider signature from this 
framework. We stick to have a high-scale universal Aq set to zero. For all the points sgn{fi) 
has been taken to be positive and the Higgs mass parameters has been set equal to the third 
generation scalar masses rn'jj^ = rnj^^ = mo^^. The values of tan /3 has been chosen such that 
it satisfies the experimental constraint branching ratio for the b — > 57 [19] which at the 3a 
level is 

2.77 X 10"^ < Br{b S7) < 4.33 x 10"^ (3) 

Parameters are fine-tuned in a way that it gives a correct cold dark matter relic abun- 
dance. In cases it is smaller than the WMAP data [20], which at 3cr is 

0.091 < ncDMh^ < 0.128 , (4) 

where ^Icdm is the dark matter relic density in units of the critical density and h = 0.71 ± 
0.026 is the reduced Hubble constant (namely, in units of 100 km s~^ Mpc~^). 

^For breaking through 770, we quote the resuh, when it breaks through the Pati-Salam gauge group 
G422D iSU{4)c X SU{2)l X SU{2)ii with even D-parity and assumed to break at the GUT scale itself. 
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This leaves us with the opportunity for some other dark-matter component. As we 
mentioned earher we compare our result with a mSUGRA point, which has a similar gluino 
mass. We denote the benchmark points as BPl, BP2 and BPS, while the mSUGRA point 
is denoted as MSGl. We tabulate the high-scale and low scale parameters in Table 2. BP2 
has the lightest gluino among these and BPS has the heaviest one. 

Note that for BPS, we do not need a gaugino mass non-universality to keep the lighter 
chargino, neutralinos to rise above the experimental bound. 

For renormalization group equation ROE, we use the code SuSpect v2 . 3 [21] and stick 
to two-loop RGE with radiative corrections to the gauginos and squarks. We use full one 
loop and dominant two loop corrections for the Higgs mass. The low scale value of the strong 
coupling constant has been chosen at asiMz) = 0.1172. We ensure radiative electroweak 
symmetry breaking and the electroweak symmetry breaking scale has been set at y/rrif^ rrif^ , 
the default value in the code SuSpect. We compute the cold dark matter relic density with 
the code microOmega v2.0 [22]. 

3 Collider Simulation and Results 

We now discuss the collider signatures of the benchmark points advocated in the preceding 
section. 

We first discuss the strategy for the simulation which includes the final state observables 
and the cuts employed therein. In the next subsection we discuss the numerical results 
obtained from this analysis. 

3.1 Strategy for Simulation 

The spectrum generated by SuSpect v2.S as described in the earlier section, at the benchmark 
points are fed into the event generator Pythia 6.4.16 [2S] by SLHA interface [24] for the 
simulation of pp collision with centre of mass energy 7 TeV and 14 TeV. 

We have used CTEQ5L [25] parton distribution functions, the QCD renormalization and 
factorization scales being both set at the subprocess centre-of-mass energy All possible 
SUSY processes and decay chains consistent with conserved -R-parity have been kept open. 
We have kept initial and final state radiations on. The effect of multiple interactions has been 
neglected. However, we take hadronization into account using the fragmentation functions 
inbuilt in Pythia. 
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parameter 


BPl 


BP2 


BP3 


MSGl 


tan /3 


17 


12 


23 


5 


(M3,M2,Mi) 


(140,180,180) 


(125,250,250) 


(153,153,153) 


(200,200,200) 


(mo^mo^'^) 


(160,1000) 


(160,1000) 


(160,1000) 


(70,70) 


















+ 


+ 


+ 


+ 


fi 


166 


123 


190 


285 


m-g 


389 


349 


421 


495 




1037 


1034 


1043 


457 




187 


142 


217 


324 




364 


348 


382 


501 




268 


253 


285 


427 




302 


273 


330 


447 




1005 


1012 


1003 


159 




157 


178 


138 


107 




104 


98 


98.5 


132 




223 


244 


233 


317 




225 


245 


232 


319 


X3 


178 


135 


203 


291 


X2 


110 


125 


101 


134 


m^o 

Xi 


64 


71 


56 


73 




0.08 


0.035 


0.128 


0.128 




3.04 X 10-^ 


4.3 X 10-^ 


4.27 X 10"^ 


3.09 X 10-^ 



Table 2: .■ Benchmark points BPl, BP2, BPS and MSGl(masses are in GeV). The first five 
parameters define the model, while the rest are low scale prediction. 
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As we have mentioned earlier that the spectra we have chosen for study is best discovered 
with a 6-reach final state. The final states studied here are: 

• (46) : 46 + X + Efr 

• (26 +2^)■. 26 +2^ +Ejr 

where £ stands for final state isolated electrons and or muons, Efr depicts the missing energy, 
X indicates any associated lepton or jet production. 

We will discuss these objects in details, that constitute the final state observables. The 
nomenclature assigned to the final state events in parentheses will be referred in the following 
text. 

All the charged particles with transverse momentum, px > 0.5 GeV^ that are produced 
in a collider, are detected due to strong B-field within a pseudorapidity range |?7| < 5, 
excepting for the muons where the range is |?7| < 2.5, due to the characteristics of the muon 
chamber. Experimentally, the main 'physics objects' that are reconstructed in a collider, are 
categorized as follows: 

• Isolated leptons identified from electrons and muons 

• Hadronic Jets formed after identifying isolated leptons 

• Unclustered Energy made of calorimeter clusters with pt > 0.5 GeV (ATLAS) and 
|?7| < 5, not associated to any of the above types of high-E'r objects (jets or isolated 
leptons). 

Below we discuss the 'physics objects' described above in details. 

• Isolated leptons {iso i): 

Isolated leptons are identified as electrons and muons with pt > 10 GeV and \ri\ <2.5. 
An isolated lepton should have lepton-lepton separation AR^^ >0.2, lepton-jet separation 
(jets with Et > 20 GeV) AR^j > 0.4, the energy deposit Yl due to Xosn-Et hadron 
activity around a lepton within Ai? < 0.2 of the lepton axis should be < 10 GeV, where 
Ai? = \/ /\rf + A0^ is the separation in pseudo rapidity and azimuthal angle plane. The 
smearing functions of isolated electrons, photons and muons are described below. 

^This is specifically for ATLAS, while for CMS, pt > I GeV is used. 
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• Jets {jet): 

Although we do not use any exphcit jet veto in our specific final states excepting for the 
b-tagging we would like to discuss the jets as they are one of the important ingredients of 
the final state observables in a collider. Jets are formed with all the final state particles after 
removing the isolated leptons from the list with PYCELL, an inbuilt cluster routine in Pythia. 
The detector is assumed to stretch within the pseudorapidity range |?7| from -5 to +5 and is 
segmented in 100 pseudorapidity (r/) bins and 64 azimuthal (0) bins. The minimum Et of 
each cell is considered as 0.5 GeV, while the minimum Et for a cell to act as a jet initiator 
is taken as 2 GeV. All the partons within Ai?=0.4 from the jet initiator cell is considered 
for the jet formation and the minimum Ylparton -^t"'^* for a collected cell to be considered as 
a jet is taken to be 20 GeV. We have used the smearing function and parameters for jets 
that are used in PYCELL in Pythia. 



We identify partonic h jets by simple 6-tagging algorithm with efficiency of = 0.5 for > 
40 GeV and |r/| < 2.5 [26]. 

• Unclustered Objects {Unc.O): 

As has been mentioned earlier, all the other final state particles, which are not isolated 
leptons and separated from jets by Ai? >0.4 are considered as unclustered objects. This 
clearly means all the particles (electron/photon/muon) with 0.5 < E^ < lOGeV and |?7| < 5 
(for muon-like track |?7| < 2.5) and jets with 0.5 < E^ < 20GeV and |?7| < 5, which are 
detected at the detector, are considered as unclustered objects. 

Once we have identified the 'physics objects' as described above, we sum vectorially the x 
and y components of the momenta separately for all visible objects to form visible transverse 
momentum {pT)vis, 



where, Y^Px = Y.iPx)iso i + Y.{Vx)jet + Y.[Vx)unc.o and similarly for Y.Vy. We identify the 
negative of the {pT)vis as missing energy E^. 



• b-jets: 




(5) 



E^ 



(Pt) 



VIS 



(6) 



Finally the selection cuts that are used in our analysis are as follows: 
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• Missing transverse energy Efr > 100 GeV. 

• Pt^ > 20 GeV for all isolated leptons. 

• Et^-^^' > 40 GeV and \r]jet\ < 2.5 

Interestingly the channels we study, specifically the ones associated with Ab-s, do not have 
any dominant SM background at all. Although, the 46 production at LHC is significantly 
high, it is expected to get down with the high missing energy associated with the final states. 
For the case with 2b, however, tt production gives the most serious background. We generate 
SM events in Pythia for the same final states with same cuts. Our cuts and event selections 
have been motivated by [27]. 

Studies in similar direction have also been discussed in references [28]. 

3.2 Numerical Results 

Here we tabulate the event rates obtained for different channels at different benchmark points 
advocated above. 



Model Points 


Total 


99 


iiil 


BPl 


141.1 


50.5 


45.3 


BP2 


315.6 


95 


164.5 


BPS 


100.6 


33 


23.6 


MSGl 


93.7 


12 


2.5 



Table 3: Total supersymmetric particle production cross-sections (in pb) as well as some 
leading contributions for each of the benchmark points. Ecm= H TeV. 

In Table 3 and Table 4, respectively for Ecm=14: TeV and 7 TeV, we note the total 
2—7-2 supersymmetric production cross-section, as well as the dominant processes, namely 
gluino and ti pair production. We can see BP2 has the highest production cross-section and 
then comes BPl, BP3 and MSGl respectively. This can be easily attributed to the mass 
hierarchy of the SUSY mass spectra chosen for these benchmark points. For BP2, the stop 
pair production cross-section is even more than the gluino production at 14 TeV. Cross- 
sections for all the benchmark points decrease significantly at 7 TeV. Most interestingly, 
gluino production cross-section decreases significantly at 7 TeV for all the benchmark points. 
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Model Points 


Total 


99 




BPl 


19.9 


3.9 


6.0 


BP2 


48.1 


8.1 


26.0 


BPS 


16.0 


2.2 


2.6 


MSGl 


12.4 


0.6 


0.2 



Table 4: Total supersymmetric particle production cross-sections (in pb) as well as some 
leading contributions for each of the benchmark points. Ecm= 7 TeV. 



Benchmark Points 




o'm 


0'2h2l 


BPl 


1.50 


0.15 


0.43 


BP2 


1.23 


0.16 


0.46 


BPS 


1.17 


0.17 


0.33 


MSGl 


0.04 


0.01 


0.40 


ttbar 


0.00 


0.00 


0.26 



Table 5: Event-rates (pb) in multichannels at the chosen benchmark points for Ecm= 14 
TeV. CTEQ5L pdfset was used. Factorization and Renormalization scale has been set to fip = 
= V^j subprocess centre of mass energy. 

This is understandable as the gluon flux decreases at 7 TeV to a large extent which reduces 
the gluino production. Hence, for BPl, BP2 and BP3 the stop production cross-section is 
more than the gluinos at 7 TeV. 

We also point out that for BPl, gluino dominantly decays to sbottom-bottom, while 
sbottom dominantly decays to bottom neutralino. The stopl in this case, decays dominantly 
to bottom-chargino. This remains the case for BP2 and BP3 as well. However, for the 
mSUGRA point, gluino decay is equally divided among all the squarks. Hence, we can easily 
see that MSGl lags far behind in channels advocated here, except for the 2b +21 + E/r 
channel. This is clear from Table 5 and 6, where we note the cross-sections for individual 
channels at Ecm= 14 TeV and Ecm= 7 TeV respectively. This is also due to the production 
cross-sections as noted earlier. 

From Table 5, we see that the number of events in channel Ah + X + E^ is as large as 
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Benchmark Points 




0"4M 




BPl 


0.115 


0.013 


0.029 


BP2 


0.080 


0.0129 


0.0423 


BPS 


0.0765 


0.0099 


0.0234 


MSGl 


0.0017 


0.00062 


0.0311 


ttbar 


0.0000 


0.0000 


0.0420 



Table 6: Event-rates (pb) in multichannels at the chosen benchmark points for Ecm= '7 TeV. 
CTEQ5L pdfset was used. Factorization and Renormalization scale has been set to fip = fifi = 
y/I, subprocess centre of mass energy. 

1500 even for luminosity lfb~^ for 14 TeV for BPl. The number of 46 + i + events 
is also as large as 150 for the same luminosity at this point. For BP2 and BP3, these 
decrease, but still they appear with events of same order of magnitude. Contribution from ti 
background in these channels are almost zero. Hence, these are really promising channels to 
discover SUSY at the LHC The channel 2b +2i + E^ although has a contribution from ti 
background, the significance {S/ y/B, where S and B signify signal and background events) 
for BPl is as large as 26 at Ifb"^. For BP2 and BP3 the significance is almost of the same 
order, which is sufficiently over the discovery limit. We also note that, as expected, MSGl 
although has large number of events to be seen at this energy, but significantly smaller than 
the non-universal benchmark points particularly in channels 46 + X + E^ and 46 + £ + E^. 

At 7 TeV, although the number of events decrease significantly, still number of 46 +X +E/r 
events for BPl is as large as 115 at 1/6^^. All the features discussed for 14 TeV, remains 
the same for 7 TeV as well. This indeed warrants our claim for an early SUSY discovery for 
the model and benchmark points advocated in our analysis. 

We also show distributions in different channels for (i) Missing Energy, (ii) Sum over 
Pt of 6s, leptons with Missing energy. The last one is really close to the definition of so- 
called Effective Mass. We show these distributions at 14 TeV. Missing energy distribution 
mostly indicate the LSP mass of the underlying benchmark points. This distribution peaks 
around 2m -n. Hence, the flatter distribution is for MSGl among the benchmarks. On the 
other hand, the effective mass distribution picks out the mass of the dominant production 
process. The peak of the distribution is around 2m§. Here as well, MSGl shows the flattest 
distribution among others. While these distributions show the robustness of our analysis. 
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it is also true that this can't be used really to distinguish these models, which have masses 
quite near to each other. However, that is not our intention really in this analysis. 
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Figure 1: Missing energy distribution in 4b, 4bl and 262/ events from left to right and top to 
bottom before putting cut. CTEQ5L pdfset was used. Factorization and Renormalization scale 
has been set to fip = fJ'R = V^, sub-process centre of mass energy. 



4 Summary and Conclusions 

In this work, we have proposed a scenario in which discovery of supersymmetry may be 
possible at the early runs of the LHC. Our scenario includes a light gluino and light 3rd 



14 




200 400 600 800 1000 1200 1400 1600 1800 2000 

Sum Of pt of b and I and MET(GeV) 



0.12 r 



0.1 - 



o 
Z 



c O.Oi 



LU 
O 

o 



0.04 



Z 0.02 





600 800 1000 1200 1400 1600 1800 
Sum of pt of b and I and MET(GeV) 



BPi 

BP2 

BP3 

MSG1 

ttbar 



,1^ jy^JHjri^ ^^' ' I .T ^ ,- „ 



200 400 600 800 1000 1200 1400 160 

Sum of pt of b and I and MET{GeV) 



1800 2000 
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from left to right and top to bottom before putting cut. CTEQ5L pdfset was used. Factorization 
and Renormalization scale has been set to iip = i^R = V^, sub-process centre of mass energy. 
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family squarks, whereas the squarks for the first two famihes are very heavy. We work within 
the gravity mediated SUSY breaking scenario, but go beyond mSUGRA using non-universal 
boundary conditions for the gaugino masses, as well as for the sfermion masses. The gluinos 
are light enough to be pair produced copiously at the LHC. Their decays to sbottom-botom 
and stop-top gives rise to bottom quark rich final states with or without charged leptons 
together with large missing energy associated with the lightest (LSP) neutralino. The most 
promising final states are Ah+E/r, Ah+l+E/rand 2b+2i+E^. We have chosen three benchmark 
points and also one mSUGRA point with light gluino mass for comparison. With the cuts 
appropriate for the LHC, we find that the signal cross section for the 4b + channel is well 
above 1, 000 fb, whereas for the other two channels, it is well above 100 fb. For comparison, 
the signal for the mSUGRA point (MSGl) is much smaller. The dominant background with 
the used cuts are from ti production which is significant for the 2b + 2£ + Ejr channel, but 
negligible for the 46 -|- Ejr, Ab + i + E/r channels. Thus as the LHC accumulates data this year 
to few hundred pb~^ (or lfb~^), the prospect of discovering SUSY in our proposed scenario 
is excellent., and we urge the ATLAS and CMS Collaborations to search for SUSY in these 
final states. 

NOTE ADDED: During the completion of writing of this manuscript an article [29] in 
similar direction appeared in the arXiV. These authors also emphasize similar final states as 
in our work. They do not specify any SUSY breaking model, but rather use three different 
branching ratio for the gluino decays to choose their benchmark points. In our work, we use 
gravity mediated SUSY breaking framework, and generate the SUSY spectra from the high 
scale with non-universal boundary conditions, and satisfy the experimental constraints on 
the SUSY particle masses , dark matter and other low energy processes. Our results for the 
cross sections for the final states are in good agreement with their results. 
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